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Abstract. Acoustic propagation in a moving fluid provides a conceptually clean 
and powerful analogy for understanding black hole physics. As a teaching tool, 
the analogy is useful for introducing students to both General Relativity and fluid 
mechanics. As a research tool, the analogy helps clarify what aspects of the physics 
are kinematics and what aspects are dynamics. In particular, Hawking radiation is 
a purely kinematical effect, whereas black hole entropy is intrinsically dynamical. 
Finally, I discuss the fact that with present technology acoustic Hawking radiation 
is almost experimentally testable. 
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Physics and Cosmology. (Springer-Verlag) . 

1 Developing the analogy 

To ask how sound waves propagate in a moving fluid is a surprisingly sub- 
tle question that rapidly introduces one to the full power and complexity of 
curved-space Lorentzian differential geometry ||l|,||J|,^,||j6l0] . A sound wave 
propagating in a flowing fluid shares many of the properties of a minimally 
coupled massless scalar field propagating in a non-flat (3-1-1 )-dimensional 
Lorentzian geometry. This partial isomorphism is the basis of a very useful 
analogy whereby parts of General Relativity can be identified with parts of 
non-relativistic fiuid mechanics. Kinematic aspects of GR, such as the exis- 
tence of event horizons, carry over to fluid mechanics (event horizons map 
into the boundaries of regions of supersonic flow). Dynamic aspects of GR 
(the Einstein equations) do not carry over. The analogy is not an identity, 
nevertheless enough features are shared in common to make the model very 
useful, and rather entertaining. (Since this is a summer school, I will be very 
pedagogical and will set out a number of exercises as we work through the 
details.) 

1.1 Ingredients 

The basic idea is to consider a non-relativistic, irrotational, barotropic fluid. 
The fluid should be irrotational since in this case the velocity is completely 
specifled by a scalar fleld, (which does not have to be single- valued): 

Vxv^O; V ^ Vi). (1) 

Thus there is hope that the sound waves, which we shall soon see are merely 
linearized fluctuations in the velocity field, can also be described by a scalar 



2 Matt Visser 



field. There is of course nothing physically wrong with non-scalar sound (vec- 
tor sound?) but the equations then become so unmanageable as to be com- 
pletely unwieldy. In the meantime, as long as whatever vorticity happens to 
be present is confined to thin vortex tubes, the present analysis is perfectly 
capable of handling everything outside the vortex core. The fluid should be 
barotropic: The pressure should be a function of the density only. This re- 
quirement makes sure that the pressure forces do not generate vorticity — it 
guarantees that an initially irrotational fluid will remain irrotational. An ad- 
ditional simplifying assumption, for the purposes of this talk, is that the vis- 
cosity is zero (inviscid fluid) . This is merely a technical simplification, and the 
complications attendant on introducing viscosity into the system are already 
understood (Viscosity acts as an explicit breaking term for the acoustic 
Lorentz invariance, and acoustic Lorentz invariance becomes an approximate 
symmetry that improves in the low- frequency low-wavenumber limit.) The 
relevant dynamical equations are: 

(1) the continuity equation 

|^ + V.(H = 0; (2) 

(2) the (zero-viscosity) Euler equation 

p(^^ + iv- V)^^) = -Vp; (3) 

(3) some barotropic equation of state 

p = p{p). (4) 

Exercise: (Easy) Check that this is a closed system of equations. Check that 
if the vorticity is initially zero it will remain so. 



1.2 Manipulations 

Now pick some arbitrary but exact solution [po(i, x), po(t, x), ipoit, x)] of the 
equations of motion. Treat this exact solution as a background field and ask 
how linearized fluctuations around this background behave. Write 



p{t,x) = pQ{t,x) + e pi{t,x) ^ , (5) 

p{t,x) = po{t,x) + epi{t,x) ^ , (6) 

tp{t,x)=Mt,x)+eiJi{t,x) + ---. (7) 

The equations of motion for these linearized fluctuations are 

^ + V- (piVVo+PoV^i) = 0, (8) 

Po (^^+V^o-V^i) (9) 

Pi = cIp,. (10) 
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Here we define = dp/ dp, the standard definition for the speed of sound. 

Exercise: Derive these equations. (The second equation, which conies from 
the Euler equation, is a Uttle tricky.) 

These three first-order partial differential equations can be assembled into 
one second-order partial differential equation which carries exactly the same 
information: 



d ( 2 (cli'x 



(11) 



This is a second-order partial differential equation for •\\}\ with variable co- 
efficients that depend only on the background field around which we are 
linearizing. Once a solution Tpi{t,x) has been obtained, the pressure fluctua- 
tions pi{t,x) and density fluctuations pi follow directly from the linearized 
Euler equation and the linearized equation of state. 
Exercise: (Easy) Check this. 

To turn this into a form suitable for obtaining a spacetime interpretation, 
introduce four-dimensional coordinates via the usual identification 



{t,x). 



Now introduce a 4 x 4 matrix 



1 



PoCs 



-1 



(12) 



(13) 



(Greek indices run from 0-3, while Roman indices run from 1-3.) Define 

g = [det ignr' . (14) 

Then the rather formidable- looking second-order partial differential equation 
for ^pl can be very simply written as 



1 d 



-g ^x^^ 



-g 9 



d 



^1=0. 



(15) 



Exercise: (Straightforward) Check this. Calculate g. Show that ( p^ ) above 
is identical to (0) . 



1.3 Interpretation 



Once you have reduced the equations of motion to the form (|T^), the last 
step is trivial: Just observe that this equation is exactly that of a minimally 
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coupled massless scalar field propagating in a spacetinie with inverse metric 
g'^'^{t,x). In fact, the differential operator appearing in this equation is just 
the d'Alambertian of the inverse metric g'^'^{t,x). 
Exercise: Check that the matrix 



is the inverse of 5^''. Check that the signature of this matrix is (— , +, +, +). 

Thus we have demonstrated that the propagation of sound is governed 
by an acoustic metric — Sfiuit, x). This acoustic metric describes a (3 + 1)- 
dimensional Lorcntzian (pseudo-Riemannian) geometry. The metric depends 
algebraically on the density, velocity of flow, and local speed of sound in 
the fluid. This is the essential difference between this acoustic Lorcntzian 
geometry and GR: The acoustic metric is governed by the fluid equations of 
motion (continuity, Euler's equation, and the barotropic equation of state) 
which constrain the background geometry, and the Einstein equations of GR 
are not useful in this context. You can certainly calculate the Ricci tensor 
and Einstein tensor for this acoustic metric, but there is no justification for 
asking these quantities to satisfy any particular constraint. 

Note that although the underlying physics (fluid mechanics) is completely 
non-relativistic, sharply separating the notions of space and time, the fluctu- 
ations (sound waves) nevertheless couple to a spacetime metric that places 
space and time in a unified framework. 

Exercise: (Some tricky points.) Copy/extend the standard definitions of 
black hole, event horizon, apparent horizon, and surface gravity into this 
context. Check your ideas against the discussion below, and the more detailed 
formulation in 

2 Examples 
2.1 Nozzle 

A particularly simple example of a non-trivial acoustic geometry is provided 
by laminar fluid flow through a nozzle : As the nozzle narrows the fluid 
speeds up. If the nozzle is sufficiently narrow the fluid velocity will exceed the 
local speed of sound. (Doing this experimentally while maintaining laminar 
flow is a rather difficult proposition.) If the fluid velocity exceeds the speed of 
sound, then sound waves cannot escape back out of the region of supersonic 
flow. Thus a region of supersonic flow shares many of the properties normally 
associated with a black hole (more properly, this region shares many of the 
properties of the ergosphere of a black hole). 




(«o)i 



gf,u[t,x} = — 



(16) 



(wo)i 
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2.2 Vortex 

A draining bathtub, with water swirling around the drain, is a useful model 
for emphasizing the difference between an event horizon and an ergosphere Q . 
As one moves inwards towards the drain, two interesting things happen: First, 
the magnitude of the fluid velocity exceeds the speed of sound, and second 
(somewhat nearer the drain) the radial component of the fluid velocity exceeds 
the speed of sound. 

The region in which the fluid velocity exceeds the speed of sound defines 
the ergoregion. (ft is impossible to stand still in the ergoregion without pro- 
ducing a sonic boom.) Once the radial component of the velocity exceeds the 
speed of sound, then acoustic disturbances cannot escape from the region 
around the drain, and this defines the event horizon. 

ff the motion is perfectly radial, (no swirling) the two notions agree. It is 
only if there is a swirling motion near the drain that the two notions need 
to be separated. This is the analog for fluid dynamics of the GR behaviour 
of the metric near a rotating black hole: For the Kerr metric the dragging 
of inertial frames implies that the region in which one cannot remain at rest 
with respect to asymptotic infinity [the ergoregion] is not the same as the 
region from which you cannot escape to asymptotic infinity. 



2.3 Supersonic Cavitation 

It is experimentally very easy to set up a situation in which air bubbles in 
water are induced to collapse at supersonic speeds. (Speeds of up to Mach 
4 are quite common.) Supersonic bubble collapse provides an example of an 
apparent horizon (not an event horizon). It's an apparent horizon because 
simply by waiting for bubble collapse to stop, and the re-expansion phase 
to start, you can always be guaranteed of getting a sound signal back out 
to spatial infinity — thus there cannot be a true event horizon (absolute 
horizon) in the system. 

Furthermore, this apparent horizon can exist for only a very short time 
during each collapse cycle: By construction, the apparent horizon lasts for 
less than one sound-crossing time. 

Experiments of this type are normally set up to investigate the phe- 
nomenon of sonoluminescence. Before anyone gets too carried away, let me 
state explicitly that the visible light emitted in sonoluminescence is not Hawk- 
ing radiation associated with this apparent horizon: (f) If anything, you 
should expect phonons, not photons (2) The "Hawking temperature" 
estimated from the acceleration of the bubble wall, while somewhat larger 
than the most naive estimates based on Unruh's analysis Q], is still far too 
small to be relevant for sonoluminescence Ql; (3) The fact that the apparent 
horizon lasts for less than one sound-crossing time renders Hawking's calcu- 
lation moot. See below for more discussion of acoustic Hawking radiation. 
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2.4 (Conformal) Schwarzschild flow 

Can we find a fluid flow tliat exactly mimics the Schwarzschild geometry? No, 
but we can get reasonably close: We can find a fluid flow that has an acoustic 
metric that is conformal to that of Schwarzschild spacetime. Start by writing 
the Schwarzschild geometry in Painleve-GuUstrand form 



ds^ = - (^1 - dt^±^^^ dr dt+dr^+r^ {dd^ + sin^ 9 dc^^) . (17) 

In this coordinate system the Schwarzschild geometry has been written in 
such a way that space is flat, though spacetime is curved. 

Exercise: Find the coordinate transformation needed to go from any of 
the more usual representations of Schwarzschild spacetime to this one. 

Exercise: Demonstrate that for any spherically symmetric geometry (not 
necessarily static, though you may want to consider that special case first) it 
is always possible to find a coordinate system such that space is flat. (So all 
the spacetime curvature can be forced into the gu and gu components of the 
metric.) 

Exercise: Take the general acoustic metric. Pick Cs a position-independent 
constant, v = ^ 1GM / r f, and p oc r"'^/^. Check that this fluid flow satisfies 
the equation of continuity Find the equation of state. Find the background 
pressure distribution needed to satisfy the Euler equation. Demonstrate that 
for this choice of fluid flow 

(^//t^)acoustic OC r ^ (^/^i^) Schwarzschild- (-^^) 

Finally, show that this is the best that can be done [p|. 



3 Surface gravity 

In the same way that one can define a surface gravity for a black hole in 
General Relativity, it is also possible to set up a notion of surface gravity for 
an acoustic black hole. Unruh showed that under certain conditions the 
surface gravity is related to the normal derivative of the fluid velocity as it 
crosses the event horizon, and is then equal to the physical acceleration of 
the fluid as it crosses the event horizon. 

dv 

gn = Cs = "fluid- (19) 
on 

Unfortunately, this result is limited to the case where (1) the speed of sound 
is position independent, and (2) the fluid crosses the event horizon perpendic- 
ularly (which means the event horizon must be identical to the ergoregion). 
The general result, derived in ||], is 
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So a position dependent speed of sound, or a non-trivial ergoregion, greatly 
complicates life. 

Exercise: There are a large number of technical incantations required to 
justify these formulae. Consider a null geodesic that just skims the acoustic 
event horizon, and parameterize it by non-relativistic Newtonian time. Show 
that this parameter is not an affine parameter for the null geodesic, and show 
that the surface gravity measures the extent to which Newtonian time fails 
to be an affine parameter. For more details see M. 



4 Acoustic Hawking radiation 



With the build up we have seen so far, the discussion of acoustic Hawking 
radiation pP,P,p^ is almost anticlimactic, and can be relegated to a series 
of exercises. 

Exercise: Read the original paper demonstrating the existence of Hawk- 
ing radiation (the one in Nature compare with jTsj). Check that this 
derivation does not need or use the Einstein equations. (The key feature of 
this original derivation is that the black hole is quasi-static: there should be 
an apparent horizon that lasts for a long time compared to the light-crossing 
time for the black hole.) 

Exercise: Compare this with some of the subsequent rederivations of the 
Hawking radiation effect. For example, using analytic continuation to Eu- 
clidean signature. Later derivations are technically slicker (and more compu- 
tationally efficient) but often obscure the underlying physics. 

Exercise: Verify that an acoustic black hole will emit a quasi-thermal 
phonon spectrum with temperature 



kT, 



H 



27rc, 



2^ 



{Cs - VA_) 



(21) 



Near the horizon, the spectrum is almost exactly thermal. As the phonons 
move away from the event horizon they are to some extent back-scattered by 
the acoustic metric. Exactly the same phenomenon occurs in GR and is the 
origin of the famous grey-body factors — even for a Schwarzschild black hole 
the emission spectrum is not exactly Planckian. 
Exercise: Put in some numbers. Verify that 



Th = 1.2 X 10"^ K mm 



1 km/sec 



1 d{cs - VA_) 
Cs dn 



(22) 



Thus for supersonic flow of water through a 1 mm nozzle, Th ~ 10~^ K. If 
this number was just a little bit better, we could reasonably hope to build lab- 
oratory experiments to verify this acoustic Hawking effect. Temperatures of 
10~^ Kelvin are not by themselves completely out of reach, though you would 
certainly not be using water as the working fluid. The real issue is that of 
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detecting a thermal phonon spectrum at this temperature, while maintaining 
laminar supersonic flow for the working fluid. 

Exercise: Verify that the existence of Hawking radiation is a purely kine- 
matic effect. Hawking radiation will occur in any Lorentzian geometry that 
contains an event horizon, independent of what the dynamical equations un- 
derlying the geometry are. See for instance [^6|jl|,|l|,|^,|o|j2^,||,|3|j2|j2|,||. 

5 Horizon entropy 

In contrast to Hawking radiation, which is a purely kinematic effect, the 
notion of black hole entropy is intimately tied up with the dynamics of the 
geometry. In fact, if the dynamics of the geometry is governed by a Lagrangian 
that depends only on the metric (and thus implicitly on the Riemann tensor 
of the metric) then there is a rather general formula for black hole entropy 

^ = i / e,, V^g d'x. (23) 



el J 5Ra 



Here the integral runs over the two-dimensional event horizon, and e de- 
notes the two-dimensional Levi-Civita symbol defined on the event horizon. 
In particular, suppose (for simplicity) that we have 



'C=^i? ( l+l^a„i?" ) , (24) 
then for a spherically symmetric black hole 

(25) 




This example is enough to drive home a key point: Entropy equals area (plus 
corrections) if and only if the Lagrangian is Einstein-Hilbert (plus correc- 
tions) |2^,^ . So calculations of black hole entropy are implicitly calculations 
of the Lagrangian governing the geometry Q . 

There are some interesting quirks of history here: Historically Bekenstein's 
notion of black hole entropy came first, and Hawking radiation was discovered 
as a side effect of trying to make the notion of black hole entropy consistent 
with ordinary thermodynamics. But now we see that Hawking radiation is a 
much more primitive concept, one that is more fundamental than the black 
hole entropy it helped explain. In fact Hawking radiation makes perfectly 
good sense even in situations in which the notion of black hole entropy is 
entirely meaningless j?]. 

This also has implications for string theory Q : We have known since the 
mid 1980's that the low-energy (sub-Planckian) limit of essentially any string 
theory is a theory of curved spacetime with dynamical equations derived 
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from an action that is the Einstein-Hilbert action (plus corrections). Thus in 
any phenomenologically interesting string theory black holes must have an 
entropy that is proportional to the area (plus corrections). In a certain sense, 
complicated state-counting calculations in the underlying string theory can 
be viewed as consistency checks that verify that the low-energy dynamics is 
what you thought it was. There are suspicions, though not a complete proof, 
that it might be able to formalize this statement by rephrasing the state- 
counting calculations directly in terms of the low-energy degrees of freedom. 
See for instance the Horowitz-Polchinski "Correspondence Principle" 
or Carlip's and Solodukhin's analysis in terms of a the central charge 
of an appropriate Virasoro algebra attached to the event horizon. 



6 Discussion 



The acoustic model for Lorentzian spacetime is a very good toy model for 
forcing you to think long an hard about fundamental issues in GR (and fluid 
mechanics, and even String Theory). It forces you to sharply separate those 
aspects of black hole physics that are purely kinematical from those parts that 
are intrinsically dynamical. It forces you to think about the universality of 
Lorentzian geometry: Even completely non-relativistic fluid dynamics has a 
Lorentzian spacetime hiding inside it. It allows you to formulate in a coherent 
manner possible approaches to the breakdown of Lorentz invariance (though 
I have not said anything about this topic in this talk). 

A key result that I would like the reader to appreciate is this: Hawking 
radiation is kinematics; Black hole entropy is dynamics. 

Finally, an observation: It is often quite remarkable how much really deep 
and fundamental physics can be found hiding in quite unexpected places. 
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